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Static Shape Control of Spacecraft 
Structural Elements Using 
Piezoelectric Actuators 


Abstract— This paper discusses piezoelectric actuators which are used in high precision 
structures in order to maintain static shape such as spacecraft antenna reflectors, sensor 
base structures, telescope etc. for its optimal performance, which is very challenging task. 
Surface errors are often introduced by pre-stress and thermal distortions due to 
temperature differences. In aerospace industry metals are replaced by composite 
materials for their improved structural properties, but when these materials subjected to 
large variation of temperature in space also yield small distortions because of Coefficient 
of Thermal expansion (CTE), which is different for constitutive materials present in its 
layer, even after using thermal insulations. So for these kinds of applications, smart 
structures involving attached piezoelectric actuators have been proposed earlier [1, 2]. 

A representative model is taken to study the application of piezoelectric materials in static 
deformation control. Linear piezoceramic materials are considered in this study. 
The behavior of piezoelectric material alone or in combination with other materials was 
analyzed by using Finite Element Method (FEM) through MSC® Marc Mentat software. 
Based on these simulation results, it was possible to better understanding the behaviour of 
piezoelectric material and thus to predict its optimum location, voltage range to obtain a 
better performance. This paper includes analysis results from other papers in order to 
validate our procedure and capability of MSC® Marc Mentat software; it also includes 
the contribution of different piezoelectric strain coefficients on the displacement in 
piezoelectric analysis. The developed simulation methodology and modeling can be 
applied to other types of linear piezoelectric materials. 


Index Terms—static shape control (SSC), piezoelectric material, finite element 


method, PZT, MSC® Marc Mentat. 


I. INTRODUCTION 


Static shape control (SSC) is the manipulation of the 
structure by means of actuation, to guide the structure to 
conform to a desired shape [1]. Static shape control is 
mainly applied in flexible beams and large space structures. 
In stealth technology, Shape Control may be implemented 
to reduce protuberance or to decrease radar cross-sections 
of aircraft. Other Shape Control applications in space 
structures include an inflatable deployable space antenna, a 
tendon control system for flexible structures and a space 
antenna reflector. Some of these space _ structure 
applications must consider thermal effects in their Shape 
Control procedures.The application of static shape control 
was investigated in this paper particularly for aluminium 
and composite plate and antenna reflectors configuration 
using piezoelectric actuators. In first section of this paper, 
in order to validate our results comparisons are made with 
some literature and theoretical results. In second section 
mainly focuses on the shape control of structures when 
subjected to thermal distortion using piezoelectric 
actuators. Also includes some results to give optimum 
location of actuators, best suitable material and voltage 
range for high actuation or shape control. 
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II. VALIDATION 


A. Piezoelectric Single layer subjected to voltage 


A Piezoelectric Single layer of size 76x26x1 mm is 
subjected to voltage change or potential difference at top 
and bottom surfaces and the FEA results obtained are 
compared with Analytical results. This FEA plate model 
consists of 76x26x1 numbers of Hex 8 elements (Element 
type 163 in MSC Marc Mentat). Material properties used 
for this analysis is PZT-5H is given in TABLE A. 
Boundary conditions used for this validation study are, 
100V is applied at top face (Z=0.001m), OV is applied at 
bottom face of element (Z=0 m) and mechanical constrains 
are on y-axis, U=w=0 and on x -axis, v=w=0. 


Hand Calculations are made using below formulae 


d33 Vh 591x107'*x100x0.001 
WwW = SS 





= 5.91x10° 
h 0.001 elles 
— dg3VW  —274x107**x100x0.026 | Sse 
a 0.001 Se ase 
d3,VL —274x10~!*x100x0.076 4 
u= — = _ = -2.11x10-°m 


h 0.001 


TABLE I COMPARISON BETWEEN THE FINITE ELEMENT RESULTS AND 
ANALYTICAL SOLUTION 
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The doming of THUNDER Piezoelectric actuator when 
eee janet Me | ae subjected to different voltage change or potential 
Marc Mentat) Error differences across the thickness are obtained for all the four 
groups and is summarized into a table below. 
direction) 


TABLE IV SUMMARY OF RESULTS FOR ALL GROUPS 


Voltage Applied | Displacement in thickness(Z) direction in mm 
(V) 
0.03643 | 0.03506 0.1133 0.1059 


v (width direction) 6.99x107 -7.124x10" 
u (length direction) | -2.11x10° 2.082x10° 


It 1s observed from Table I that the perfect match oe ee eae 





4 0.07286 | 0.07011 0.2267 0.2118 
between analytical and FE analysis results exist in_ the 0.09108 1 0.08764 0.2833 0.2648 
thickness direction or across poling direction, but slight 0.10930 | 0.10520 0.3400 0.3177 


variation of FE results in other two directions mainly due to 


boundary conditions From TABLE IV, it is concluded that higher actuation 
is possible for Group 3, where thickness of piezoelectric 


B. Validation of High Displacement Piezoelectric Actuator backing material 1s less and in plane dimensions large 
(THUNDER) Finite Elements Models (compared to other Groups) for same PZT-5A material and 


input voltage value. 





NASA’s Langley Research Centre [2], have conducted 
both experiments and analyses on THUNDER piezoelectric 
actuator to predict doming during the manufacturing 
process and subsequent straining due to an applied input 
voltage. Analysis results are obtained by ANSYS and 
NASTRAN. This study is conducted in order to compare 
our results with Experimental and analysis results. Based 
upon thickness and length of base metal (steel), they 
divided into four groups, for these group both analyses and 
experiments are conducted separately. Tables below give 
material properties and geometric model configurations. 


Group 2 


Group 1 


Displacement, mm 
72.8 


Displacement, mr 





Voltage,V Voltage, 


TABLE If COMPOSITION OF TEST SPECIMENS 


TABLE I PROPERTIES OF ALL SPECIMENS (ALL GROUPS) 


: / , 
Layers Material inc’ Modulus _ /|Poisson's |Coefficients e icn | Experimental | & 15641 Y _— 
10° m of Ratio |of Thermal 106 01 ve 106-1 — 
Elasticity Expansion 50602 06.0 
10° rC 50 150 250 350 ie a | 2 m mw 


Voltage,V 
Aluminum 2: p 254 68.95 0.33 24 
Layer) Figure 2. Comparison of Experimental and Different FEA results 


| 2 | LaRC-SIE | 254 | 4 | 045 | 46 

Also From fig.2, graphs show the comparison for 

different groups and it is clear from these graphs that 

coe | eae ee | ae | oe | ae | MSC Marc Mentat software package produces linear 
results and these results closely associated with 

Time: 2.000e+00 = eae _ experimental results when compared to ANSYS and 
NASTRAN. 


Group 3 Group 4 





Displacernent, mm 
~~ we we 


oa 
=f 
Displacement, mm 









i 


177e-04 
860e-04 


542e-04 


C. Thermal Distortion Analysis 


dh In this type of analysis, irrespective of structures 
material properties, when it is subjected to change in 
temperature, it will deform depending on the Coefficient of 
Thermal Expansion and distance from fixed end. This 
— validation is to show and calculate the thermal distortions 
— fy due to change in temperature for simple space antenna 
Br cael — ‘ reflectors of different materials, which are needs to be very 
precise for optimal performance. 
I) Space antenna reflector of 0.6m diameter and 


F/D=0.4 


906e-04 
588e—-04 
271e-04 


527e-05 





349e-05 


weno F EF FN NN ®N 


| 
“NI 





Figure 1. Displacement of Fourth Group specimen for 300Volt 
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MEL Softecwre 


Space antenna reflector of diameter 0.6m and F/D (F- 


Focal Length and D- diameter of reflector) ratio 0.4 is | 
subjected to 10°C change in temperature and various se 
displacements are calculated and tabulated below for . 
different materials like Aluminium, Kapton, and Composite > 
etc. | 
Hand Calculations ‘ 
a) For Aluminium cote 
Coefficient of Expansion of Aluminium a = 23x10° °C ee 
Displacement X = a x AT x L - 


=23x10°x10x0.3=6.9x10° m 





Displacement Y=ax els L : Figure 3. Displacement in Y (thickness) direction for Space KAPTON 
= 23x10° x10x0.06 = 1.38x10~ m antenna reflector 
Displacement Z = ax AT x L Il. ANALYSIS 


— -6 — <5 
= 23x10" x10x0.3 = 6.9x10" m This section consists of shape control of cantilever beam 


and Space antenna reflector, using piezoelectric actuators 


b) For KAPTON when they are subjected to thermal distortion. two types of 





Coefficient of Expansion of KAPTON a= 20x10° PC analyses are carried out, thermal distortion analysis due to 
temperature change and piezoelectric analysis due to 
Displacement X = 0.x AT x L application of voltage or potential difference to 
= 20x10° x10x0.3 = 6.0x10° m piezoelectric actuators. For cantilever beam or plate these 
analysis are carried out, in order to find out best location, 
Displacement Y = a x AT x L permissible thickness of actuators, voltage or potential 
= 20x10° x10x0.06 = 1.20x10° m difference for piezoelectric actuators for better actuation or 
shape control. First section is to find out best location and 
Displacement Z = 0x AT x L thickness and in second section, all those concepts are 
= 20x10° x10x0.3 = 6.0x10° m applied in shape control of space antenna reflector 
subjected to thermal distortion due to change in 
c) For Composite Material temperature. 
TABLE V COMPOSITE MATERIAL PROPERTIES A. Cantilever beam 
The analyses are carried out on cantilever beam to 
Gist oe 7 find best location, permissible thickness of actuators, and 
voltage or potential difference for piezoelectric actuators. 
These parameters for best actuations are found out on 
simple cantilever beam made up aluminium and composite 
Density = 1717kg/m° materials, when it is subjected to a unit degree change in 
Displacement X = 0),;XATXL temperature across it thickness. 
— 5.2x10° x10x0.3 = 1.56x10° m Six different cases are made and analyzed by 
Displacement Y = 033xATXL changing position of actuators, and thickness piezoelectric 
— 21.6x10° x10x0.06 = 1.29x10° m patch. Even though the depoling voltage is depends on 
Displacement Z = o.xATxL thickness, type of piezoelectric material used, but in these 
— 5.2x10° x10x0.3 = 1.56x10° m analyses up to 400V is applied across the electrode which 


need not to be applicable in experiments because high 
TABLE VI SUMMARIES OF THEORETICAL RESULTS .; , 
electric field (generally depoling will starts over 


Mitedals CTE | Length | Height Chan : Disp x | Disp y | Disp z 1000V/mm electric field). 
°C (m) (m) a (m) (mn) an) 
Aluminium io aia 6.900" | a) Aluminium Cantilever Beam attached with piezopatch 


subjected to temperature change 
/KAPTON | ~ | 6.00e" | O0e° | 1.20e° | 6.00e 4 ; ’ 
| 9.06 | 10 | 6.006% “ssa The aluminium cantilever beam of length 100mm, width 


rene 5 Tone [we Pre Tanase 10mm and thickness Imm is subjected to 1°C change in 
temperature across its thickness. Because of temperature 
difference the cantilever beam will deflect. To regain the 
original shape of cantilever beam, piezoelectric actuators 
are used. Piezoelectric patch of dimension 20mm (length), 
10mm (width) and by varying thickness with constant 
voltage value percent regain in tip deflection is measured. 





FE Analysis results are same as above table 
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Placing Unimorph and bimorph of PZT patches of 1mm,or 
0.75mm, or 0.5mm or 0O.25mm thick at a position from 
fixed end are either x=0, 10, and 20 mm. PZT -5H material 
properties used in this analysis. 


It shows from table that, composite beam requires very 
less DC voltage to regain its tip deflection compared to 
aluminium cantilever beam, because thermal distortion 1s 
less in composite when compared to isotropic aluminium 


material. These above concepts and conclusions are applied 
to reflector in next section. 


B. FEA for reflector surfaces 


In this section two types of analysis are done on 
structural elements used in space industry i.e. Space 
antenna reflector. Thermal load is considered as a main 
reason for distortion of space structures and this distortion 
may be of very less in magnitude. In order to regain shape 
of precise structures piezoelectric actuators are used. So in 
this section, both thermal distortion analysis first and 


TABLE VII VOLTAGE REQUIRED TO REGAIN 100% TIP 


DEFLECTION 


Aanainiaecantleves 
Position From fixed end [Position From fixed end 
Beam 
(mm) (mm) 
PZT \x=0 |x=10 |x=20 | x=o |x=10 |x=20 
thickness 


a 
ae 
—_ 


Voltage 


required 
for 100% 





* Theoretical values and these values may not be applicable in experiments. 


From above summary, the maximum actuation is 
possible when the piezoelectric patch is placed near the 
root i.e. where the maximum strain density present and 
thickness of piezoelectric patch should be very less. 
Bimorph require very less voltage value compared to 
Unimorph, because bimorph itself add more stiffness to 
structure compared to Unimorph and there are two 
oppositely polarized Piezopatches mounted on either side 
surface of the cantilever beam. 


b) Composite beam attached with piezopatch subjected to 
temperature 


The composite cantilever beam of length 100mm, width 
10mm and thickness Imm is subjected to 1°C change in 
temperature across its thickness. Because of temperature 
difference the cantilever beam will deflect. To regain the 
original shape of cantilever beam, piezoelectric actuators 
are used. Piezoelectric patch of dimension 20mm (length), 
10mm (width) and by varying thickness with constant 
voltage value percent regain in tip deflection is measured. 
Placing Unimorph and bimorph of PZT patches of Imm,or 
0.75mm, or 0.5mm or 0.25mm thick at location from fixed 
end are either x=0, 10, and 20 mm. 

Material Properties used in this analysis are 

Elastic moduli (GPa): E);=E..=17.2 E33=6.9 

Shear moduli (GPa): Gj2=1.7 G,3=G23=2.76 

Poisson’s ratio: y;2=0.14 y13=y23=0.4 

Density p=1717 kg/m” 

Thermal expansion coefficient (um m! °C'): o);=0.=5.2 
033=21.6 


TABLE VIII VOLTAGE REQUIRED TO REGAIN 100% TIP 


DEFLECTION 


Position From fixed 
end (mm) 


Position From fixed 
end (mm) 


Composite Plate 


Voltage 


Required 


noe X=10 | X=20 | X=0 | X=10 | X=20 
thickness 


piezoelectric analysis second are carried out in sequence 
for same model. 





Figure 4. FE models of Reflector surfaces 


1) Aluminium reflector surface 

Following Finite elements of MSC Marc Mentat have 
been used in the FE modeling. The details of these 
elements are given below: 


ELASTIC 
8 noded, HEX, Element type 7, ACTIVE Degrees of 
Freedom, Ux, Uy and Uz (#3 translations) 


PIEZOELECTRIC 

8 noded, HEX, Element type 163, ACTIVE Degrees of 
Freedom, Ux, Uy and Uz and v (#3 translations, | electric 
potential) 


a) Analysis of 0.6m dia., aluminium reflector with 12 
numbers of PZT-5H patches 


FE investigation has been carried out for 0.6m dia. 
(F/D=0.4) Aluminium (0.5mm thick) parabolic shell with 
12 numbers of piezo patches, PZT-5H using MSC Marc 
Mentat software, to predict the deflection pattern of 
reflector. 

All three linear translational DOFs of the reflector 
have been restrained in the center of the reflector. The 
patches are affixed on the convex side of the surface along 
the main X-Z axis, 90 degree apart with 3 patches on 
each side of the centrally held point for the reflector. Both 


yoy, Le 
oo 
05mm [135 [100] 81} 82} 68 | 130 


thermal distortion analysis for 10°C and _ piezoelectric 
analysis for 200V are carried out and results are plotted 
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TABLE Ix DISPLACEMENT IN THICKNESS DIRECTION FOR 


ALUMINIUM REFLECTOR 


Analysis type Displacement 


Thermal Distortion 
Piezoelectric actuation 


Max. Reflector tip displacement FEA is 5.596x10°’m 
(0.56 microns). Less displacement is attributed to higher 
stiffness of the small size of the reflector. 

For this case not much actuation is possible because 
thermal distortion is in the order of 10’s microns and 
piezoelectric actuation is in 0.1’s microns. 





b) Analysis of 0.6m dia. Aluminium reflector with 4 
numbers of PZT-5H patches 


FE investigation has been carried out for 0.6m dia. 
(F/D=0.4) Aluminium (0.5mm thick) parabolic shell with 4 
numbers of piezopatches, PZT-5H using MSC Marc 
Mentat software, to predict the deflection pattern of 
reflector. 


All three linear translational DOFs of the reflector have 
been restrained in the center of the reflector. The patches 
are affixed on the convex side of the surface along the main 
X-Z axis, 90 degree apart with 1 patch on each side of the 
centrally held point for the reflector. Both thermal 
distortion analysis for 10°C and piezoelectric analysis for 
200V are carried out and results are plotted. 


MSC Software 








FE model on MSC Marc Mentat for 0.6m dia. Aluminium 
reflector with 4 Unimorph 


Figure 5. 


MSL A Sottonare 











Figure 6. Thermal distortion of 0.6 m dia. Al. Reflector with 4 piezo 
patches subjected 200 V 
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Figure 7. Compensated Thermal distortion of 0.6 m diameter Aluminium 
Reflector by 4 piezo patches subjected 200V 


Thermal distortion is compensated by one order with 
only 4 Piezopatches attached near fixed end. So this 
concept 1s applied to regain shape of composite material 
reflector in the next section. 


2) Composite reflector surface 

FE investigation has been carried out for 0.6m dia. 
(F/D=0.4) composite (0.08mm thick) parabolic shell with 4 
numbers of piezo patches, PZT-5H using MSC Marc 
Mentat software, to predict the deflection pattern of 
reflector. 

All three linear translational DOFs of the reflector 
have been restrained in the center of the reflector. The 
patches are affixed on the convex side of the surface along 
the main X-Z axis, 90 degree apart with | patch on each 
side of the centrally held point for the reflector. Both 
thermal distortion analysis for 10°C and _ piezoelectric 
analysis for 200V are carried out and results are plotted. 

Four piezoelectric layers each 1s having thickness of 0.2 
mm are considered in this analysis. 


Material Properties used in this analysis are 

Elastic moduli (GPa): E,;=E.=17.2 E33=6.9 

Shear moduli (GPa): Gj2=1.7 G)3=G23=2.76 

Poisson’s ratio: Y;2=0.14 y13=y23=0.4 

Density p=1717 kg/m” 

Thermal expansion coefficient (um m'! °C"): O11 =O9=5.2 
033=21.6 


The details of these elements are given below: 

ELASTIC 

8 noded, HEX composite, Element type 149, 
ACTIVE Degrees of Freedom, Ux, Uy and Uz (#3 


translations) 


PIEZOELECTRIC 

8 noded, HEX, Element type 163, ACTIVE Degrees of 
Freedom, Ux, Uy and Uz and v (#3 translations, | electric 
potential) 
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Figure 8. Thermal distortion of 0.6 m dia. composite Reflector with 4 
piezo patches 
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Figure 9. Piezoelectric actuation of 0.6 m dia. Composite Reflector with 
4 piezo patches subjected 200V 


TABLE X DISPLACEMENT IN THICKNESS DIRECTION FOR COMPOSITE 


REFLECTOR 


Analysis type Displacement 


thermal Distortion 11.59 microns 
piezoelectric actuation | 38.62 microns 





Max. Reflector tip displacement by thermal distortion is 
1.159E-05m (11.59 microns) and by piezoelectric actuation 
1s 3.862E-5 m (38.62 microns) 

From this case, the shape of reflector surface can be 
regained by only 4 Piezopatches with less than 200V DC. 


IV. CONCLUSIONS 


Studied the Electromechanical behavior of Piezoelectric 
(PZT) material in the actuation, Validation study is carried 
out to establish a procedure to analyze structural models 
attached with piezoelectric actuators and compared results 
with theoretical and experimental results. It is found that 
the result shows close correlation with the test results 
compared to the other software. 

Thermal deformation analysis using Finite element 
analysis of different structural elements is carried out. 
Cantilever beam with different configurations are studied 
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for tip displacement control. It 1s demonstrated through 
Finite element analysis that the small displacement can be 
controlled by using piezoelectric materials. 

An attempt is made to use piezoelectric material for 
reflector. Reflector of diameter 0.6m is considered, with 
Aluminium as well as composite. Different PZT locations 
are studied. Since the thermal deformation in composite 
reflector is very less compared to the aluminium reflector, 
deformation compensation shows better results in 
composite reflector. 
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TABLE A PZT -5H PIEZOELECTRIC MATERIAL 
PROPERTIES 


Material: PZT-5H (Anisotropic) 
Density=7500kg/m 

Elastic coefficient matrix(C):6x6 N/m 
C),=Cx=1.26x 10" Cy2=Co:=7.95 x 10°” 


Dielectric or permittivity matrix (k):3x3 Farad/meter 
€11=€2=1.505x10"; 633=1.301x10™ 


iezoelectric Strain coefficient properties (d):3x6 C/N 


2 1=d39=274x 10 


15=d24=74 1x10 


33=091x10 
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